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Abstract 



The Airy gas model of the edge electron gas is used to construct an 
exchange-energy functional which is an alternative to those obtained in 
the local density and generalized gradient approximations. Test calcula- 
tions for rare gas atoms, molecules, solids and surfaces show that the Airy 
gas functional performs better than the local density approximation in all 
cases and better than the generalized gradient approximation for solids 
and surfaces. 
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Since the pioneering papers on density functional theory (DFT) |],|2] there has been a 
constant search for exchange-correlation functionals of chemical accuracy. This includes 
the works on the generalized gradient approximation (GGA) |3|-[7| which are dedicated 
efforts to construct local functionals for inhomogeneous systems ranging from atoms to 
solids based on the uniform electron gas, i.e., the local density approximation (LDA), 
and density gradient corrections, as well as the development of a number gradient 
level, semiempirical functionals f5|-|TT||. The GGA functionals have had a considerable 
impact upon the fields of quantum chemistry and solid state physics because they 
reduce the LDA overbinding and generally improve the calculated properties, relative 
to experiments, of molecules [p|,p!2|-p^| and bulk solids |T5|-pO||. However, they perform 



less well for the bulk properties of late transition metals and semiconductors [piH 23|l, and 
the underestimate of the exchange energies of surfaces as well as the overestimate of 
the dissociation energies of the multiply bonded molecules p|,[T3|| indicate the necessity 
to go beyond the gradient level approximations and develop functionals that depend 
upon other inhomogeneity parameters, e.g., higher derivatives of the charge density or 
the Kohn-Sham kinetic energy density. One step in this direction is the meta generalized 
gradient approximation (meta-GGA) of Perdew, Kurth, Zupan, and Blaha |12| which 
proves highly promising for both finite and extended systems . 

In the present work we introduce and apply a new gradient level exchange energy 
functional based on the concept of the edge electron gas ||24j| . Besides the formal in- 
terest in the development of density based, orbital independent functionals there are 
several reasons why in applications of DFT the focus is on the approximate, local 
exchange-correlation schemes. Within the Kohn-Sham approach to DFT the Kohn- 
Sham exchange energy may be determined exactly and as demostrated recently [p5 -27 



so may the corresponding local exchange potential. However, the exact Kohn-Sham 
exchange formalism is non-local and orbital-based, i.e. both the exchange energy and 
potential are highly complicated non-local functionals of the Kohn-Sham orbitals. In 
consequence, the application of exact exchange is computationally demanding. Further- 
more, when exchange is treated exactly the error cancellations between the exchange 
and correlation energies on which all approximate schemes depend are lost owing to the 
poor description of correlation effects and, as a result, the total energies worsen fl5| . |26 |. 
For these reasons the exact Kohn-Sham exchange energy has only been used in practice 
in connection with semi-empirical, hybrid approximations [ [2~8|j . 

The concept of the edge electron gas was put forward by Kohn and Mattsson |23| as 
an appropriate basis for the treatment of systems with edge surface outside of which all 
Kohn-Sham orbitals decay exponentially. Its simplest realization, the Airy gas model, is 
based on the linear potential approximation and may serve as the starting point for the 
construction of functionals which are alternative to the GGA. The Airy gas model has 
recently been used to construct an explicit kinetic energy functional for inhomogeneous 
systems |29| which for atoms and surfaces has the accuracy of functionals based on a 
second order gradient expansion. 

Here we have taken the exchange energy of the Airy gas model derived by Kohn and 
Mattsson |24]] and cast it in a form amenable to a simple, accurate parametrization. The 
procedure may be viewed as local mapping of the real system described by its density 
and scaled gradient onto the Airy gas model and represents one possible solution to 
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the joining of the interior to the edge regions. The parametrized functional which we 
refer to as the local Airy gas (LAG) functional is tested in calculations of the exchange 
energies of rare gas atoms and of metallic surfaces within the jellium model where the 



exact results are known [p0|| . In addition, we apply the LAG exchange functional in 
conjunction with the LDA for correlation in calculations of the molecular binding 
energies and bulk properties of solids. 

The present LAG exchange functional has a number of advantages over previous 
GGA functionals: i) it explicitly includes the properties of the edge region where much 
interesting physics occurs, ii) its accuracy may be systematically improved by including 
higher order expansions of the effective potential of the model system, and iii) the re- 
sulting exchange-energy functional is as simple and well-defined as that of the standard 
LDA. i.e., it has no adjustable parameters. 

The starting point for the Airy gas exchange energy functional is the potential 

/ \ f oo for z < —L , v 

VeffW ~ | Fz for _ L < z < ^ , (1) 

which is linear in z, independent of x and y, and has a hard wall at — L far from the 
electronic edge at z — 0. The slope of the effective potential F = dv e ff/dz leads to a 
characteristic length scale 

and the electron and exchange-energy densities are then given by 

n(z)=r 3 n(() } (3) 

and 

= -^" 4 e x ((), (4) 
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A contour plot of the universal function g(s, s') may be found in Ref. p4[ . The exchange 
energy (H) may be written in the form 



e x (z) = e L x DA (z)F x [s(z)], (7) 
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where e J x jDA {z) is the exchange energy density of the uniform electron gas. The en- 
hancement function 
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is the unique function (F x ((), s(Q) of the scaled gradient 
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plotted in Fig. |l|. For comparison we also present results obtained by the GGA of 
Perdew, Burke, and Ernzerhof (PBE) as defined in Ref. |J, and the second order 
gradient expansion (GEA) ||32|| . It follows from the figure that the exchange density 
(0) in the low gradient limit of the Airy gas model reduces to e^ DA (z) as it should. In 
the large gradient limit e x (C) ~ — , Ref. [p4H , and similar to the case of the kinetic 



energy density we use the properties of the Airy gas to find the following explicit 
asymptotic expression 



e x [n(z)} 
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in terms of the density and its derivatives. 

The density of the exchange energy per electron of the Airy gas is plotted as a 
function of the distance z from the electronic edge in Fig. 0. It is seen that the large 



gradient expression (|TD) is accurate for z/l > —1.4 corresponding to s > 0.5. It is also 
seen that neither the LDA nor the PBE GGA || leads to the correct behaviour near 
and beyond the electronic edge at z = 0. 

The scaled gradient is conserved when going from the real electron gas to the Airy 
gas model [25| and therefore the enhancement function F x (s) parametrized, for instance, 
in a modified Becke form | 33| 



^LAGi 



1 + 7S 



a\8 '■ 



which includes the proper LDA limit, may be used to obtain the exchange energy density 
of the real electron gas from the local, scaled gradient s[n(^)]. For a = 2.626712, (3 = 
0.041106,7 = 0.092070, and 5 = 0.657946 we find that the local deviation between the 
exact result (f|) and the parametrized form (pT|) integrated over the range < s < 20 
is less than 0.3%. We note that the present parametrization, being an overall fit, does 
not reduce to the GEA [|l^,0] in the low gradient limit. In contrast to the case of the 
kinetic energy |29j we have not been able to find an explicit, analytical expression for 
the exchange energy for small s values, and to establish the behaviour numerically has 
not been attempted because the s — > limit is reached only at z — » — oo as seen in Fig. 
^|. The exact behaviour of the LAG exchange functional at s — > is therefore not know 
at present. 

In the following we report the results of applying the LAG exchange functional to 
four test systems: i) rare gas atoms, ii) diatomic molecules, iii) jellium surfaces, and iv) 
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solids. In all cases the total energy is calculated using self-consistent LDA densities. For 
molecules and solids the LAG exchange energy is combined with the LDA correlation 
energy |31[, since correlation effects has not been worked out in the Airy gas model. 
The motivation of this combination is given in terms of the enhancement function over 
the local exchange energy 0, defined as F xc (s) = e xc [n]/e LDA (n), where e xc [n] denotes 
the exchange-correlation energy density. 

Most of the currently applied approximate density functionals are based on error can- 
cellations between the exchange and correlation energies [|1^,^]. For physically inter- 
esting densities this cancellation leads to F xc (s) with negligible slope up to s ~ 1. Plots 
of the enhancement function over the local exchange energy for gradient level and meta- 
GGA approximations can be found in Refs. |pjl3l. In the present LAG exchange plus 



LDA correlation scheme this function becames F^ AG {s) = F^ AG {s) + e^ DA (n) / e x DA (n ) , 
where e^ DA {n) is the correlation energy density of the uniform electron gas. Thus, the 
F£ AG (s) is determined only by the LAG enhancement function (|TT|), which, for s < 1, is 
a slowly increasing function of s. Therefore, we expect the present exchange-correlation 
scheme to preserve the excellent cancellation properties of the LDA and PBE GGA, 
and, at the same time, to bring the calculated properties in closer agreement with 
experiment than conventional LDA. 

For the rare gas atoms included in Table | the GEA, PBE, and LAG functionals 
yield exchange energies which are, on the average, 6.4 %, 8.5 %, and 1.8 %, respectively, 
larger than those obtained in the LDA. The PBE values are in very good agreement 
with the exact Kohn-Sham results |T3| , |34]| , which are given relative to the LDA energies 
in the last column of the table. The LAG approximation represents only a minor 
improvement relative to the LDA total atomic exchange energies. 

The effect of the gradient correction to the LDA atomization energies for a few se- 
lected diatomic molecules is shown in Table [Ij which also includes the relative difference 
between the LDA results and experimental data |T3. Here, the LDA charge densities 



for the molecules have been generated using the full charge density (FCD) technique in 
conjunction with the exact muffin-tin orbital method (EMTO) j35"H37f. It is seen that 
the LAG approximation (i.e. LAG exchange and LDA correlation energy) and PBE 
GGA have comparable accuracy: Both functionals reduce the LDA overbinding, and 
yield atomization energies which are, on the average, 16.8 % (PBE) and 16.2 % (LAG) 
smaller than the LDA values. 

In Fig. |3] we compare four exchange functionals applied to the jellium model of 
metallic surfaces [[3(J . The fact that for a given r s - value the exchange energies become 
increasingly negative in the order LDA, LAG, GEA, and PBE is a simple consequence 
of the enhancement functions shown in Fig. [l| and in agreement with the observation 
that the GGA significantly underestimate surface exchange energies |U. We note that 
LAG approximation represents an improvement over both the LDA and PBE and vary 
less with r s than either of the other two approximations. 

As a final test of the LAG approximation we have calculated the atomic volumes 
and bulk moduli of several metals and semiconductors in their observed low tempera- 
ture crystal structures by means of the FCD-EMTO method p5|-p7[]. The results for 
the equilibrium atomic radii are plotted in Fig. |j. For some selected metals and semi- 
conductors, for which accurate LDA, PBE, and meta-GGA results have been published 
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13fl , the atomic radii and bulk moduli are presented in Table [II. The comparison 



of our LDA |3l| atomic radii for the transition metal series with those obtained by 



the full-potential linear muffin-tin orbital [|HJ and linear augmented plane wave [|2T| 
methods using the same LDA gives mean deviations of 0.33 %, 0.43 %, and 0.49 % for 
the 3d, 4d, and 5d series, respectively. For Li and Na the present LDA results agree 
within 0.07 % with the full-potential values from Ref. [15]. We therefore expect that 
the results of the present LAG and PBE calculations shown in Fig. ^ will deviate less 
than 0.5 % from full-potential calculations. The mean deviations between the present 
atomic radii and bulk moduli listed in Table |TTT| and those of Ref. [I"3| obtained using 
the linear augmented plane wave method are 0.20 % and 3.28 % for the LDA and 0.27 
% and 3.26 % for the PBE functionals. 

The LDA atomic radii shown in Fig. [| deviate, on average, by 2.26 % from the 
experimental values [j38,39l, while those calculated in the LAG model and the PBE 
deviate by 0.83% and 0.91%, respectively. Among the energy functionals considered in 
Table |T| the LAG is found to give the lowest mean deviations for both atomic radii 
and bulk moduli. We note that for these solids the LAG approximation achieves the 



accuracy of the recently developed meta-GGA [12.13 



We have used the Airy gas model of the edge electron gas that is equivalent to the 
linear potential approximation to develop an exchange energy functional which may 
serve as an alternative to the functionals based on the generalized gradient appoxima- 
tion, e.g., PBE GGA. Test calculations for finite and extended systems show that the 
LAG approximation is more accurate than the local density approximation in all cases. 
While the LAG results for atoms are very close to the LDA results and, hence inferior to 
the PBE GGA results, its accuracy for the atomization energies of diatomic molecules 
is similar to that of the PBE GGA. In bulk systems the LAG results are, on average, 
closer to the experimental values than those obtained in the PBE GGA. These results 
are very satisfactory in view of the fact that the LAG exchange functional is derived 
solely form the properties of the Airy gas, and, hence, with no a priory assumptions 
concerning the exchange enhancement factor. In this sense it is truely ab initio but for 
the correlation effects which needs to be worked out in the Airy gas model. 
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TABLES 



TABLE I. The effect of GEA |32|, PBE §], and LAG gradient corrections (in percentage) 
on the LDA atomic exchange energies. All functionals are evaluated from the self-consistent 
LDA [|3l] Kohn-Sham densities. KS denotes the relative difference of the exact and LDA 
exchange energies from Ref. 

ATOM GEA PBE LAG KS_ 

Hi 133 Id!) 42 16\0 

Ne 6.7 9.4 1.9 9.7 

Ar 5.1 7.7 1.4 8.3 

Kr 3.5 5.5 0.9 5.9 

Xe 2.9 4.7 0.7 5.0 



TABLE II. The effect of PBE |>) and LAG gradient corrections (in percentage) on the 
LDA atomization energies for diatomic molecules. Both functionals are evaluated from the 
self-consistent LDA pi Kohn-Sham densities generated by the FCD-EMTO pHJl. Expt. 



denotes the relative difference of the experimental and LDA atomization energies from Ref. 



MOLECULE PBE LAG Expt. 

L12" -22.6 -17.8 2~X~ 

Be 2 -25.9 -32.8 -76.6 

CO -9.8 -8.6 -13.3 

N 2 -12.2 -10.7 -14.5 

NO -14.2 -11.2 -23.0 

2 -16.1 -16.0 -31.1 
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TABLE III. Theoretical equilibrium atomic radii (in Bohr) and bulk moduli (in GPa) 
for some selected solids. The present calculations have been performed for crystallographic 
a phases using the FCD-EMTO method p5|-|37[|. The results obtained by the meta-GGA of 



Perdew, Kurth, Zupan and Blaha (PKZB) and the experimental values are from Ref. [13] 
The mean absolute values of the relative errors are shown in parenthesis. 





Slda 


SPBE 


Slag 


SPKZB 


SiExpt. 


Blda 




Blag 


BpA'ZB 


^Expt. 


Na 


3.769 


3.916 


3.927 


4.019 


3.936 


8.2 


7.6 


7.3 


7.0 


6.9 


Al 


2.947 


2.989 


2.977 


2.966 


2.991 


81.2 


75.2 


76.2 


90.5 


77.3 


Fe 


2.565 


2.645 


2.604 


2.627 


2.667 


253 


178 


209 


198 


172 


Cu 


2.602 


2.684 


2.656 


2.656 


2.658 


193 


137 


157 


154 


138 


Pd 


2.846 


2.916 


2.883 


2.888 


2.873 


235 


184 


203 


181 


181 


W 


2.929 


2.977 


2.953 


2.946 


2.940 


312 


292 


299 


311 


310 


Pt 


2.888 


2.943 


2.916 


2.908 


2.892 


304 


244 


268 


267 


283 


Au 


2.998 


3.081 


3.043 


3.041 


2.997 


194 


134 


156 


153 


172 


Si 


3.163 


3.198 


3.189 


3.200 


3.182 


100 


92.8 


94.0 


93.6 


98.8 


Ge 


3.303 


3.384 


3.354 


3.349 


3.318 


71.6 


61.2 


64.0 


64.6 


76.8 


GaAs 


3.296 


3.375 


3.346 


3.347 


3.312 


73.0 


62.0 


72.1 


65.1 


74.8 


NaCl 


3.202 


3.346 


3.337 


3.284 


3.306 


32.9 


23.0 


21.7 


28.1 


24.5 




(1.48 %) 


(1.28 %) 


(0.80%) 


(0.88 %) 




(17.2%) 


(9.2%) 


(9.1%) 


(9.3 %) 
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FIGURES 



FIG. 1. The exchange energy enhancement function (^) and the parametrized form (|TT|), 
the latter indicated by dots, of the local Airy gas (LAG) compared to those of the local density 
approximation (LDA), the the generalized gradient approximation (PBE) as defined in Ref. 
||, and the second order gradient expansion (GEA) 



FIG. 2. The exchange energy per electron, —e x (z)/n(z), of the Airy gas, obtained from 
the enhancement factors shown in Fig. |l|, plotted as a function of the distance from the 
electronic edge and compared to the exact result obtained from (|j) and (^|) and the explicit 
large gradient limit (|T0|). Energy in units of — (e 2 /2)/ _1 and distance in units of I defined in 
(|2|). The scaled gradient s is also shown. 



FIG. 3. The exchange energy of the LDA g|, GEA Q, PBE g, and LAG approxima- 
tions obtained from the self-consistent Kohn-Sham densities relative to the exact exchange 
energy [ 40 1 in the jellium surface model for a range of density parameters r s = (3/4-7rn) 1//3 . 



FIG. 4. Relative deviations of the calculated and experimental (3£| equilibrium atomic 
radii for the alkali, alkaline earth, 3d, 4d, and 5d transition metals using LDA PBE || 
and LAG energy functionals. The numbers in parenthesis are the average deviations. 
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